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Abstract 

A theory for wave mechanical systems with local inversion and translation symmetries is developed 
employing the two-dimensional solution space of the stationary Schrodinger equation. The local sym¬ 
metries of the potential are encoded into corresponding local basis vectors in terms of symmetry-induced 
two-point invariant currents which map the basis amplitudes between symmetry-related points. A uni¬ 
versal wavefunction structure in locally symmetric potentials is revealed, independently of the physical 
boundary conditions, by using special local bases which are adapted to the existing local symmetries. 
The local symmetry bases enable efficient computation of spatially resolved wave amplitudes in systems 
with arbitrary combinations of local inversion and translation symmetries. The approach opens the 
perspective of a flexible analysis and control of wave localization in structurally complex systems. 

PACS numbers: 03.65.-w, 01.55.+b 
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I. INTRODUCTION 


Symmetries play an essential role for the structure and predictions of modern physical theories 
by their generic relation to constants of motion. In classical dynamics, continuous symmetries 
lead to the conservation of associated currents following from Noether’s theorem [1] which has 
subsequently been generalized in various ways [2-5]. In a quantum description, the relation 
between symmetry and conservation laws is extended to discrete symmetries [6, 7] by the com¬ 
mutation of the corresponding operators with the Hamiltonian, thus yielding a connection to 
the possible form of stationary eigenstates of a system. In particular, states of definite parity in 
inversion-symmetric systems and conserved quasimomenta in structures with discrete transla¬ 
tion invariance (to be referred to as parity and Bloch theorems, respectively) are central to the 
treatment and understanding of a large class of phenomena in, e.g., atoms or crystals. 

The significance of symmetries is perhaps most appreciated when they are broken [8], either 
explicitly at the level of the equations of motion or spontaneously by the system state itself 
[9] . Symmetry breaking is thereby commonly related to emergent effective interactions [8] or to 
(ground) state properties [10], a prominent example being the origin of particle mass in the Higgs 
mechanism [11]. Regarding spatial transformations, the Hamiltonian of a composite system 
may obey a symmetry only in a subpart of configuration space, in which case the symmetry is 
broken globally. This restricted occurrence of a spatial symmetry constitutes a kind of symmetry 
breaking which is in fact unavoidable due to the finite size of any actual system. In reality 
any symmetry of the effective potential describing a system is indeed restricted to some finite 
spatial region, while multiple symmetries may occur domain-wise (see Fig. 1 for an illustration 
of a composite system described by a potential with different symmetries in different domains). 
Such ‘local’ spatial symmetries [12] may occur inherently in complex systems such as large 
molecules [13-15], in quasicrystals [16-18], or even in partially disordered matter [19, 20]. They 
are also often present by design in, e.g., multilayered photonic devices [21-23], semiconductor 
superlattices [24], acoustic waveguides [25, 26] or magnonic systems [27]. In such artificial setups, 
broken global symmetry is often required to obtain structures suitable for specific applications. 
A special case are completely locally symmetric (CLS) setups, where the active region is covered 
exclusively by domains with local symmetries [28]. 

Despite their omnipresence, local symmetries and their consequences in wave mechanical 
systems are largely overlooked when passing from global to no symmetry: Although it is very 
common to treat a composite structure in terms of its coupled subparts, their eventual local 
symmetries are seldom encoded in the description or directly exploited in calculations. A first 
step towards an approach addressing local symmetries was taken in Ref. [29] by defining local 
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FIG. 1. Local symmetry decomposition of a potential region into inversion (17) or translation (T) 
symmetric domains T> ( ] of size Xd — Xd-\ which are further decomposed into Nd cells C^ 1 . Nd . The 
77-symmetric domain V 2 (with Nd = 2 cells) is mapped onto itself under inversion through a, while 
in the T-symmetric domain (with Nd ^ 2 cells) V 3 the first Nd — 1 cells (denoted by V 3 ) are mapped 
onto the last Nd — 1 cells (denoted by P 3 ) under translation by L. Domains without symmetry (like, 
e. g., V 1 and V 4 ) may in general also be part of the potential region. 


inversion operators and relating them to one-dimensional (ID) wave scattering via associated 
locally, i. e. domainwise, invariant quantities. These symmetry-induced invariants, which have 
the form of two-point currents, were subsequently used to classify perfect transmission states 
in terms of their spatial profile [30]. The generic spatial structure of stationary states in the 
potential regions of local inversion or translation symmetry was recently established [28] by a 
mapping relation 




Q^i'x) - Qi>*(x) 


( 1 ) 


between the wave amplitude i/j at symmetry-related points x, x, where the complex two-point 
currents 


Q = Ji ~ ^'(a^O)], (2) 

Q = Yi - il>*\x)il>{x)] (3) 


are invariant, i. e. spatially constant, within the corresponding domain of local symmetry. We 


here use the notation ^'(x) = 


!(=.\ _ dip(x) 


dx 


Their values depend on the details of the potential via 


"0 but are related to the globally invariant current J by 


\Q \ 2 - \Q? = J 2 


( 4 ) 


in Hermitian systems. 

The above mapping relation generalizes the parity and Bloch theorems to systems where 
reflection and translation symmetries, respectively, are realized only domain-wise [28]. In fact, 
Q contains information on how a global symmetry is broken: it vanishes in the case of global 
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symmetry of both the potential and the boundary conditions, in which case Eq. (1) can be written 
as the corresponding well-known eigenvalue problems (see Ref. [28]). Q becomes nonzero and 
globally constant for asymmetric boundary conditions. In the case of a local symmetry holding 
in a certain domain, Q and 0 are constant in this domain. Interestingly, 0 remains invariant 
even in the presence of (locally) symmetric complex potentials (in contrast to the usual current 
J), as verified experimentally for CLS acoustic waveguides [31]. It also proves suitable as an 
order parameter for globally time-parity-symmetric systems, as shown in Ref. [32], Focusing 
on the wavefunction mapping induced by local discrete symmetries [28], we notice that the 
applicability of Eq. (1) depends on the boundary conditions imposed on the stationary state 0, 
since the mapping coefficients diverge for J = 0. While there is always some finite transmission 
in usual ID scattering settings, the current typically vanishes for any bound eigenstate (as well as 
for scattering eigenstates of the inversion operator [29]), and therefore an equally valid symmetry 
mapping for such cases is desirable. 

In the present work we develop a formalism for locally symmetric wave mechanical systems 
incorporating the above concepts in a form which is independent of the boundary conditions 
imposed on the physical setup at hand. Formulated directly in the (two-dimensional) solution 
space of the stationary Schrodinger equation, the theory uniformly encodes local inversion and 
translation symmetries of the potential into corresponding local basis vectors. In particular, a 
construction scheme for a global basis of an arbitrary one-dimensional wave mechanical system is 
deviced, exploiting the presence of multiple local symmetries in the underlying potential. Sets of 
linearly independent solutions, whose role is central for various types of Sturm-Liouville problems 
[33, 34], are hereby utilized to construct generalized two-point currents of mixed basis functions 
which are spatially constant within domains of local symmetry. These invariants establish a 
basis function mapping between symmetry-related points in each local symmetry domain which 
incorporates the potential symmetry in the domain basis without explicit reference to the spatial 
dependence of its components. The local symmetry bases are then matched to assemble a global 
basis in terms of arbitrary initial local basis solutions defined in a subspace (cell) of each domain 
and the associated invariant currents. Boundary conditions are finally imposed on a physical 
solution expanded in this global basis. The introduced procedure is in contrast with the usual 
approach for solving stationary Schrodinger (or Helmholtz) equation in systems with symmetries 
in restricted domains, where: (1) it is necessary to assume that within a domain the basis is 
the same as that of a system with the respective global symmetry and (2) the explicit form 
of the domain basis is used to construct the overall solution with appropriate matching. The 
introduced formalism of local symmetry bases leads to a threefold main attainment: (i) The local 
mapping relations pertaining to the generalized ID parity and Bloch theorems of Ref. [28] are 
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extended to arbitrary boundary conditions (i.e. including ones yielding zero current) in terms 
of the introduced local basis invariants, (ii) A universal structure of stationary wavefunctions in 
locally symmetric systems is revealed by their representation in local symmetry-adapted bases 
which are constructed from arbitrary initial solutions, (iii) Exploiting the diagonal form of 
the amplitude mapping matrices in the local symmetry basis (LSB) of each symmetry domain 
provides an efficient scheme for computing wavefunctions for locally symmetric systems. 

The paper is organized as follows. In Sec. II the symmetry-induced amplitude mapping 
in terms of two-function, two-point domainwise invariants is derived and the LSB, leading to 
diagonal mapping matrices, is constructed. Section III introduces the connection of LSBs of 
different domains, illustrating the generic structure of the physical solution, and demonstrates 
its efficient computation with arbitrary boundary conditions with initial input only in single unit 
cells of each domain. In Sec. IV we summarize our work and provide concluding remarks. 


II. LOCALLY INVARIANT SYMMETRY MAPPING 


Consider a potential V which is symmetric under the linear coordinate transform 

{ 2a — x (F = II) 

(5) 

x + L (F — T) 

within a domain T> C M, that is, V(x) = V(x) under an inversion (77) through a point a or 
a translation (T) by a length L (where it is understood that a T-symmetry transform applies 
to x G V° = {all but the last unit cell of V }; see Fig. 1). Any solution ip(x) to the stationary 
Schrodinger equation for this potential (setting h — m — 1), 

Hip(x ) = --^"(x) + V(x)ip{x) = Eip{x), (6) 

can be expanded in the basis of two linearly independent solutions (pi(x), 02 (:r) of Eq. (6) in V for 
a given energy eigenvalue E. Subtracting <p n (x)H<p m (x) from <p m (x) H (p n (x) , with m,n G {1,2}, 
leads to 

2 iq'mn = <t>m{x)<t>n{x) ~ K(?)4>m( x ) = 0 ( 7 ) 


within T> due to the E-symmetry of the potential. This means that the ‘mixed’ (ie., containing 
both (pi and 02 ) symmetry-induced two-point quantities 


q m n = [« 7(p m {x)(p' n {;x) - (p' rn {x)(p n {x)\ 


( 8 ) 


are spatially constant within the symmetry domain T) of the potential. In the same manner an 
alternative invariant quantity 


Qmn 


Nm(^)0n(^) - 0m(A>0n(x)] 


(9) 
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is obtained, whose translation (cr = 1) variant for m = n and x = x becomes the current 

jm = J. [C(x)0' m (>) - C(x)0 m (a;)] (10) 

corresponding to the solution (j) m . The invariants q rn n and q mfi thus have the form of mixed 
nonlocal currents, and for m = n, i. e. by replacing both <j>\(x) and 02 (%) hr Eqs. (8) and (9) 
with a single solution 0(a;), they reduce to the ‘pure’ (one-function) nonlocal currents Q and Q 
in Eq. (1). With some algebra it can be shown that the three spatial invariants q mn , q mn and j m 
are connected via the relation 

|?mn| | QVrml = jmjni (H) 

which introduces symmetry-induced constraints between the values of a single solution 0 or any 
pair of solutions 0i,02 at x and x. Before proceeding, it is worth mentioning that invariants 
analogous to the form in Eq. (8) can be derived also for more general symmetry transformations, 
as shown in the Appendix. We here restrict the presentation to the isometry transformations in 
Eq. (5), which leave the ID Schrodinger equation invariant. 


A. Symmetry mapping with mixed currents 


The mixed invariants q mfl will now be used to construct a general mapping relation between 
bases in the solution space. We first write Eq. (8) in the matrix form 


( (hi ffiA _ l_ ( M x ) 0'i(>A f> o- 0 \ / 4>\ (x) 0' 2 (x)\ 

yhl Q22/ 2l \02(a:) 4>' 2 {x)) \o -i) y0i(s) M x )) 


Left-multiplying by the inverse of the first matrix product on the right hand side and then 
transposing, we can map the solution column vector <fi = between the symmetry-related 
points x and x as 

4>(x) = Q4> (x) (13) 


via the (spatially invariant in T>) symmetry-mapping matrix 


Q 


2 i 
w 


-Q21 qi i 
\-Q22 <i\2 


1 


(14) 


where w[(j>i, 02] = 0i(^)02( a; ) — 02(ic)0 / i( a; ) lh e Wronskian of the two functions 0i, 02- If these 
are linearly independent solutions of Eq. (6) in an interval D, then w(x) ^ 0 Vx G T> [35]. Note 
that Eq. (13), with the q mn defined in Eq. (8), holds generically for arbitrary potential and for 
any coordinate transformation F. It is in the presence of 77- and T-symmetry that it becomes 
a mapping relation with constant (in T>) coefficients. For local F-symmetry with an arbitrary 
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smooth transformation F, a generalization of the bilinear mapping in Eq. (13) also exists, as 
shown in the Appendix. 

It can be shown that the Q-matrix is nnimodnlar with the determinant detQ = a distin¬ 
guishing between the case of inversion and translation symmetry. Relation (13) maps any pair 
of linearly independent solutions from x to the transformed point x in a domain via mixed 
symmetry-induced invariants. In this sense, it further generalizes the local parity and Bloch the¬ 
orems of Ref. [28] , now formulated at the level of the solution space of the Schrodinger equation 
without reference to boundary conditions imposed on physical solutions. Additionally, Eq. (13) 
indicates that 0i and 02, although linearly independent, are in general interrelated via local 
symmetry, in the sense that their values at any symmetry-related points are coupled by the 
same constant matrix Q. A decoupled pair of solutions amounts to a diagonal Q. as elaborated 
on in Sec. IIB, and is key to deriving an optimal basis for the treatment of stationary wave 
mechanical problems involving local symmetries, as discussed in Sec. III. 

For a physical solution -0 in T> (obeying the appropriate boundary conditions) which is linearly 
independent from its complex conjugate 0*, Eq. (13) reproduces the pure mapping of Eq. (1) 
(and its complex conjugate) if we choose the basis <fi to be 0i = 0, 02 = 0*- I n general, however, 
-0 and 0* are not linearly independent, as is the case for stationary bound states (which can be 
chosen real, 0 = 0*) or for stationary scattering eigenstates of the inversion operator 77 [29]. In 
these situations the current J vanishes (since it is given by the Wronskian w[ip, -0*] = 2 iJ, which 
vanishes if -0, 0* are linearly dependent solutions), and the mapping relation of Eq. (1) cannot 
be used. This limitation arises from the fact that Eq. (1) refers to a physical solution satisfying 
specific boundary conditions. The main advantage of the present approach is that symmetry- 
induced mapping relations are expressed at the level of the general two-dimensional solution 
space which is not subject to specific boundary conditions. Hence, with the basis functions 
0i,02 hr any given domain being linearly independent by assumption (so that w[0i,02] ^ 0), 
the local mapping (13) between symmetry-related points can always be exploited to construct 
a global basis as shown below in Sec. III. Expressed in this global basis a physical solution 
-0 may then have J = 0 and thus prevent the use of Eq. (1), but the local symmetries of the 
potential are already addressed in the construction of £ using the basis mappings in Eq. (13). 
In other words, with the present generalized mapping relation via arbitrary sets of independent 
domain solutions, the presence of potential symmetry is manifest in any stationary state (also 
with vanishing current) through the underlying basis: In any F-symmetric domain F, a physical 
solution = a ■ 4>(x) = ai0i(x) + a 2 02(2?) of Eq. (6) (with x G D and ai,a 2 determined by 
imposed boundary conditions) at energy E contains the symmetry information through the 
basis 0i 2 which is mapped between F-transformed points via a constant Q-matrix according to 
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Eq. (13). 


B. Local symmetry basis 


Let us now exploit the local F-symmetry in the given domain T> to arrive at a decoupled 
form of the mapping relation in Eq. (13). To this end, consider the transformation of cj) via an 
invertible constant matrix S into the special basis 


X = 



= S(f> 


(15) 


of the two linearly independent solutions x+iX- i n 77 which fulfill the mixed conditions 


ctx±(x)x' ± (x) = x!±(x)x±(x) (16) 

for any igD, that is, which has q m m = 0 (m = +,—) in V. In this basis, the invariant matrix 
performing the symmetry-induced mapping 

X (x) = Scf>(x) = SQ<f>(x) = SQS~ 1 x(x) (17) 


within T> is diagonal, 


Q x = SQS 1 = 


z .|_ 0 

0 


(18) 


since the offdiagonal elements q m fh {jn = +, —) vanish due to Eq. (16), with the eigenvalues 

trQ 


trQ / (trQ 

= Tl “ ff = 2 


± sfA 


(19) 


given by the characteristic equation z 2 — trQz + a — 0 of the original mapping matrix Q. The 
matrix S which diagonalizes Q (that is, up to a scalar factor, the inverse of the eigenvector 
matrix of Q) is given by 


S = 


7- 


-9H 


or 


v ~7+ Qil 

if q\\ 7 ^ 0 or q 2 2 7 ^ 0 , respectively, where 

1 


Q 22 -7+ 
-?22 7- 


( 20 ) 


7 ± = 2 (912 + 921 ± 1/( 9 i 2 - ’h'. r + o'C-j , 


( 21 ) 


with both matrices being equivalent if both q^j, q 22 7 ^ 0. S trivially equals the unit matrix if 
Qil — q 22 — 0 (that is, if Q is already diagonal). Recall here that the trace trQ = 2i(q 12 —q 2 l)/w = 
Z+ T Z— , the determinant a = z + z_, and thereby also the discriminant A = (trQ/2 ) 2 — a in 
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Eq. (19) remain invariant under similarity transformations in the solution space and are therefore 
real quantities, since any complex basis (p can be similarity-transformed into a real one. 

The basis x is ‘symmetry-adapted’ in V in the sense that the are eigenfunctions of the 
operator Of corresponding to the symmetry transform F acting in T>. Indeed, since Q x is 
diagonal, the LSB functions x± are n °t coupled upon their mapping, that is, each function is 
separately given by a constant factor times its image throughout the local symmetry domain T>\ 

OfX±(x) = X±(x) = z ± x±{x). (22) 

For inversion symmetry (cr = — 1 ) we always have trQ = 0 from Eq. ( 8 ) and hence A = 1, so 
that the mapping factors from Eq. (19) are 

z± = ±1 (F — 77). (23) 

Thus, x± = Xe,o is here the LSB of even and odd solutions in the domain V of a locally 77- 
symmetric potential with respect to its inversion point a. 

For translation symmetry (a = 1) we can write z + = zZ 1 = \z\e lkL , with the discriminant A 
distinguishing three cases for the z± from Eq. (19) (recall that tr Q and A are real, as explained 
above): 

(a) If A < 0, then the x± = X±k are propagating wave solutions with complex conjugate 
mapping factors 

z± = e ±ikL (F = T ) (24) 

under translation by L with kL = arctan(2-*/— A/tiQ), where |z±| J = cr = 1 is accordance with 
the conservation of each current j±. Equation (16) for the x± coincides with the condition corre¬ 
sponding to global potential symmetry for the pure mapping relation ( 1 ), that is, with vanishing 
one-function q, as shown in Ref. [28]. Therefore, k is identihed as the crystal momentum in the 
corresponding Bloch state (at energy E) for T> — M. 

(b) If A > 0, then the mapping factors can be written as real exponentials z± = e ±nL (i. e., 
in the form of Eq. (24) for imaginary k = —in) since z + Z- = cr = 1, where kL = ln(trQ/2 + 
a /A). The associated solutions x± = X,±k must now be real (up to constant phase factors) 
in order to conserve zero current under translation. Lying energetically in the gaps between 
allowed energy bands for the corresponding globally periodic system (with symmetry domain 
V — M), these solutions diverge at ±oo for k > 0, and can be involved in physically acceptable 
solutions only for setups with finite (or semi-infinite) locally symmetric domains. Enhancing 
the contribution of components of type x+k (x~k) m a physical—propagating or not—state on 
the left (right) of a boundary between local T-symmetry domains may then enable controllable 
density accumulation, that is, wave localization, around the boundary. 
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(c) If A = 0, then Eq. (19) exhibits a double root z± = 1 or —1 corresponding to k = 0 or to 
k = ±7 t/L (modulo 27 t) in Eq. (24), with associated solutions which are periodic with period L 
or 277, respectively. 


III. CONSTRUCTION OF GLOBAL BASIS AND PHYSICAL SOLUTION 

The local basis approach developed above can be used to construct a global solution basis 
£(x) for a potential with arbitrary combinations of 77- or T-symmetry domains by matching 
of the different LSB solutions at the domain interfaces, since £(x) need be continuous (and 
smooth for nondiverging potentials). The procedure enables an efficient assembly of this global 
basis £(x) which exploits the local symmetries of the system through the mapping of the LSBs 
via Q-mat rices along each symmetry domain. While the LSB mapping was formulated above 
for a single domain T >, we will now consider multiple attached domains and LSB mappings 
between consecutive cells within each domain. A labeling for domains and cells is thus introduced 
as follows: We consider a spatial decomposition of a given potential into N domains T) d = 
[xd~ i,Xd] (<7 = 1 : N = 1,2, ...,7V) which obey distinct symmetry transformations F d (that is, 
are characterized by different inversion centers a d or periods L d ). A domain T> d is further divided 
into N d cells C\= VNd of equal length (see Fig. 1). In a T-symmetry domain each cell covers a 
period L d , while a 77-symmetry domain is divided into a left and a right cell by its inversion 
point a d (so that always N d = 2 for 77-symmetry). Since the cell index has no indication of the 
domain it belongs to, each /-subscripted object is also ((^-superscripted as in C[ d \ 

We now apply the LSB mapping relation, Eqs. (17) and (18) (or, equivalently, Eq. (22)), to the 
cells of a given domain: Within a domain V d , the LSB solution x^ (obtained from an arbitrary 
basis <p (7 through Eq. (15)) is propagated from cell to cell by a diagonal mapping matrix Qx'\ 
so that 

hV/f'M) = IQfUxfV), zeef* (25) 

gives the LSB amplitude profile x\ d \ x ) di the Z-th cell of T> d in terms of the one in its first cell 
through the symmetry transform F d of the domain acting l — 1 times. 

Note here that the invariance (constancy) of the Q 1 -'^-matrix of the initial solution vector 
0(7) within a locally symmetric domain T> d allows for its diagonalization into at any pair 
of symmetry-related points in T> d - Convenient points, requiring knowledge of 0(7 in minimally 
extended regions, are the inversion point a d and the endpoints x d , x d + L of the first period for 
77- and T-symmetry, respectively. 

We proceed assuming that we have found the different LSB functions of all N domains, which 
will now be connected by matching them at the domain interfaces. As mentioned above, the 
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matching is necessary to construct a continuous (and smooth) global basis on which a physical 
solution can be represented later. The LSB in T> d will generally match a linear combination of 
the LSB solutions of the next domain V d+L at the interface x d , 

X [d \x d ) = M d+lX {d+l \x d ) } (26) 


with a matching matrix M d+ i. Note that the interface point x d may be considered to belong 
to the domain T> d or to T > d+ \, to both, or even to neither of them, depending on whether 
the symmetry transform in T> d and/or T> d+ \ applies for this boundary point or not. In any 
case, though, its left neighborhood x^ belongs to V d and its right neighborhood x^ belongs to 
V d+ 1 . Therefore, the matching conditions can be expressed in a general way using X^( x d) and 
^,(rf+ 1 )( a; +) j n (26). Special care is needed to handle the case when the potential is singular at 
x d (containing, e. g., terms proportional to 5(x — x d )), whence the matching conditions (and thus 
the matching matrix) should be adapted accordingly to the discontinuity of the wavefunction 
derivative. For a potential which allows for a continuous wavefunction derivative at x d , the 
general matching matrix reads 


M d+ i — 


w: 


d-\-\^d -|-1 


w ^ d + i w d+ + ,d 


(27) 


with shorthand notation 


= \^xf' - xflSX-z,, (28) 

following from the continuity of the functions x ( ' d \ x ) 1 M d+ iX iyd+l \x) and their first derivatives 
at x = x d . 


The aim is now to combine the above procedures—diagonal propagation of each domain’s 
LSB among its cells, Eq. (25), and matching of different LSBs at domain interfaces, Eq. (26)— 
to obtain a continuous basis for the whole potential region, denoted £(x), which consists of 
connected parts £}• d \x ) (with x G T> d ) corresponding to the different domains. Specifically, we 
start from an initial (i) desired reference domain 77*, for which we set 


£ w M = x (,) M> xev. 


(29) 


from which we shall construct the global basis ^(x) by applying matching conditions at the 
consecutive domain interfaces. Indeed, according to Eq. (26), at another domain T> d with d > i 
the corresponding part £^( x ) of the global basis will equal the LSB of that domain multiplied 
by the product of consecutive matching matrices from i to d\ 

£ (d) (a0 = | n Md '\ * (d) (- x ') = M (dl) X (d) (x), x G V d . (30) 
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However, each LSB x ^ i 11 the Z-tli cell of the corresponding domain D d can be obtained from 
the first cell through Eq. (25), with its argument back-transformed by the inverse transform F^ 1 
applied l — 1 times, 

xf\x) = [Q^^x^iF'^ix)), x G C\ d) (31) 

for all cells 1 = 1: N d . Thus, on the level of cells, the global basis can be written as a branched 
function 

£(x) = (32) 

where the (forward) basis propagation-matching matrix 

G\ di) = (33) 

first propagates x ^ from the first to the Z-th cell in D c i and then applies the matching up to 
this domain. The multidomain basis £(x) is determined by the LSB in the initial domain d = i 
in the sense of Eq. (29). The imprint of the local symmetry of the potential is manifest in £ 
through the mapping 

^ d \x) = Qf^ d \x), x£V d (34) 

within each domain via the corresponding (transformed) constant mapping matrix 

Qf = mWqWIMW}- 1 . (35) 

This reveals a universal structure of the solution space for potentials with local symmetries in 
terms of domainwise invariants. 

For / < i, the (backward) basis propagation-matching (from cell N d to l and from domain % 
to d < i) is performed by the matrix 

qM = [MM]- 1 [Q( (i) ]-^+ i (36) 

containing d — i matching matrix inversions, with the diagonal -matrices elementwise in¬ 
verted. 

If i = 1 and f = N, then £(x) constitutes a global basis for the complete potential region, on 
which the physical solution if> is expanded as 

= c ■ £(x) = Cifi(x) + c 2 6(^), (37) 

with the amplitude vector c determined by the boundary conditions imposed at x = xq,xn- As 
indicated above, the role of G^ h ' > is to propagate the LSB function x^ from the first to the 
Z-th cell in T> d and subsequently apply the domain interface matching up to this domain (see 
Eq. (32)) given an initial domain i with (^ see Eq. (29)). Thus, if the physical solution 
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has local coefficients a = af' > in an arbitrary basis (p^ (see end of Sec. II A) in the Z-th cell 
of domain Vj, then those are related to the coefficients c in the constructed global basis $, as 
a| cZ) = cG^S^i where is the matrix transforming (p^ to x ^ i n domain V^ (see Eq. (15)). 
The application of boundary conditions is thus naturally postponed until a basis of the solution 
space for the complete potential has been obtained, offering flexibility with respect to the setup at 
hand: Energy-quantizing (e.g. Dirichlet, Neumann, mixed, periodic, or exponentially decaying) 
boundary conditions determine Ci and C 2 subject to appropriate normalization, while continuous- 
spectrum (scattering) asymptotic conditions relate c li2 to propagating wave amplitudes at both 
ends (see discussion below). 

Note here the conceptual difference of the local basis approach from a conventional transfer 
matrix method where the amplitude vectors of ip are propagated in a fixed basis (usually of 
counterpropagatiug plane waves in flat potential regions): Here, instead of the physical solution, 
the basis itself is propagated in a locally symmetric setup with (repeated or inverted) unit cells 
of arbitrary potential profile. Even in cases of, e.g., intervals of finite periodic potentials, the 
corresponding basis in the interval is usually adapted from the globally periodic counterpart [36] 
with explicit spatial dependence. Here, the LSB is constructed intrinsically from an arbitrary so¬ 
lution of the first cell of the local symmetry domain. In particular, the present approach exploits 
the local symmetries by virtue of the symmetry-adapted bases which are propagated (forward or 
backward) through multiple cells by diagonal Q-matrices, thus providing an important technical 
advantage—especially in the presence of large periodic parts. 

It should be pointed out that, although the present approach is devised for potentials which 
are decomposable into multiple local symmetry domains, its application does not become invalid 
in presence of nonsymmetric domains (such as, e. g., defects in a finite periodic lattice). Specifi¬ 
cally, we can simply treat a domain T>d which is neither 77- nor T-symmetric (such as T> 1 or V 4 
in Fig. 1) as a domain with a single cell in the used notation, and use Eqs. (32) and (33) 
with / = 1 for this domain (coinciding with the cell). It is then clear that no is involved for 
this nonsymmetric domain, which makes sense since there is no local symmetry to be exploited; 
only the matching matrix M ^ remains in Eq. (33) to match this domain’s local basis Xi=i t° 
that of the previous domain. Further, since there is no symmetry to adapt the basis to in the 
domain, the basis can be chosen arbitrarily. Therefore, it is simply set equal to the initially 
computed basis, x^ — for (the first and only cell of) this nonsymmetric domain. In other 
words, the construction of the global basis £ can still be applied if some domains of the setup 
happen to be nonsymmetric, although clearly no symmetry-induced advantage can be drawn 
from these domains. 

Let us now summarize the procedure followed in the LSB approach to stationary wave systems, 
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as schematically represented by the sequence 

^(d)( C (d } ) _g l % z ^\s {d) -> X [ d) £(x) ^ (38) 

oid,L d 

and expressed as follows: 

(i) Decompose the potential into N domains T>d=i : N containing maximal regions of local 77- 
or T -symmetry, and compute an arbitrary pair of linearly independent solutions <p^ d \x ) 
(if possible analytically, or numerically with arbitrary initial conditions) to Eq. (6) only in 
the first cell C[ d ^ of each domain. 

(ii) Construct the matrix QG) from <f)( d \x) at x = a d (xd, x d + L) for local 77- (T-) symmetry, 
and diagonalize it to find its eigenvalues z^± and S',/-matrix in the basis </>0. 

(iii) Propagate and transform the first-cell LSBs Xi^ — Sdft d \z G C[ d ') within each domain 
V d by the matrices G (dt \ G ^ with a selected initial (reference) domain i and final (end) 
domains / = 1, N to obtain a global basis £(x) in the potential region. 

(iv) Impose desired boundary conditions on a physical solution 'ip(x) — c ■ £(x). 

Recall that the global basis ^ in step (iii) coincides with the local basis in the selected 
domain T>i which can be anywhere in the interaction region. Assuming that the potential can 
be, to some extent, tuned by external parameters, one could design a desired (for simplic¬ 
ity, nodeless) wave profile for x^ an d determine the corresponding—not necessarily locally 
symmetric — domain potential Vi from Eq. (6) as [37] V)(x) = E + x!±"{. x )/^X±( x )- Given the 
LSB mapping within and among domains, and in particular solutions with exponential mapping 
factors z± = e :pKi , this provides enhanced controllability of the spatial field distribution such 
as its localization in selected regions: While the coefficients cy 2 are uniquely determined in the 
case of energy-quantizing boundary conditions, in the case of scattering we can impose Cy 2 = 1 
and solve for the amplitudes a^’ > of plane waves e ±ly/ ^ x on the left (x < Xo) and on the right 
{x > Xm) of the interaction region. In other words, the ingoing amplitudes and cG that 
produce a desired domain localization at a given energy in a locally symmetric potential can be 
determined efficiently by the proposed scheme. 

IV. SUMMARY AND CONCLUSION 

Employing the two-dimensional solution space of the stationary Schrodinger equation, we 
have developed a theory for treating ID wave mechanical systems with local (i. e. domainwise) 
inversion and translation symmetries. Encoding the local symmetries of the potential into corre¬ 
sponding local basis vectors, the formalism is independent of the boundary conditions imposed 
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subsequently on particular physical solutions. The approach is based on two-function, symmetry- 
induced local invariants, which have the form of two-point currents and are spatially constant 
within each domain of local symmetry. They enable an extension of the generalized ID par¬ 
ity and Bloch theorems of Ref. [28], i.e. domainwise amplitude mapping relations, to arbitrary 
boundary conditions and thereby to states carrying zero current. More importantly, the theory 
reveals a universal structure of wavefunctions in locally symmetric potentials in terms of special 
local bases which are adapted to a given symmetry in a finite domain. The local symmetry 
bases (LSBs) are constructed from arbitrary initial solutions in only a single unit cell (one half 
of an inversion symmetry domain or one period of a translation symmetry) and mapped among 
cells by diagonal matrices. Combined with the matching of different LSBs at symmetry domain 
interfaces, this leads to an efficient computational scheme for spatially resolved wavefunctions 
in systems with arbitrary combinations of local inversion and translation symmetries. The ad¬ 
vantage of the method is especially pronounced for completely locally symmetric (CLS) systems 
with different large periodic parts. The multiplicative mapping of local basis functions within 
each locally periodic domain by exponentials then enables a natural control of wave amplitude 
distribution by tuning the potential parameters. In particular, input amplitudes in scattering 
setups can be flexibly designed to produce localization in desired spatial domains. Valid gener- 
ically for wave mechanics (e.g. acoustics, optics, or quantum mechanics), the LSB approach 
provides the perspective to use local symmetries to explain and control the amplitude response 
of structurally complex scattering or bound systems. 
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Appendix A: Local two-point invariants for general spatial transformations 

We here derive bilinear quantities analogous to the form in Eq. (8) which are spatially constant 
for general symmetry transformations of the potential of a given domain and enable a general¬ 
ization of the mapping relation in Eq. (13). Specifically, let us consider an arbitrary bijective 
coordinate transformation F : x —>• y = F(x) and a domain V mapped through F to a domain 
V , for which the potential obeys V{;x) = V ( y) with x G V and y = x G D, as shown in Fig. 2 (a). 
The aim is to construct a two-function quantity Qp(x,y ) whose total derivative with respect 
to x, in analogy to Eq. (13), vanishes under the above E-symmetry (or ‘shape invariance’) of 
the potential. With Qp being a function of y , its derivative will generally be affected by the 
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FIG. 2. Local symmetry of a potential V(x) within a domain V under a similarity transformation 
F : x —> y = F(x) (blue lines) which maps V to V with potential V(x ) for (a) a general bijective 
transformation F. (b) translation by L, y = T(x ) = x+L, and (c) inversion through a, y = IJ(x) = 2a— 
x. The locally symmetric part of the potential (colored) is generally embedded between nonsynnnetric 
parts (gray). Shaded stripes highlight the imaging of the transformed potential V{x) onto the x-axis. 


transformation F. To exploit the local E-symmetry of the potential, we therefore write the 
Schrodinger equation in the transformed coordinate y , which becomes 

Hf^f{v) = ~^D F (y) + V(y)i/> F (y) = E^ F {y) (Al) 

with the second derivative transformed to 


DfM = F"(F~'(y))Mv) + {F'(F-' (y))] 2 Mv) 


(A2) 


using the chain rule A / i/j F (y(x )) = F'(x) F (y) with x = F 1 {y) (F 1 being the inverse 
coordinate transform), where we define F'(s ) = 1 and ipF = The function i^F^y) 


denotes a solution of the modified equation arising from the deformation of the single axis of our 
ID system. Under the local shape invariance of the potential, the bilinear two-point combination 


Qf(x, y) = — i/j(x)F'(x)ip F (y) - ^\x)^ F (y) 


(A3) 


of a solution ijj(x) of Eq. (6) and a solution ^(y) of Eq. (Al) is then spatially constant in V, 


since 


2 iQ'p = il>(x)D F (y) - ^"(x)ip F (y) = 0 (A4) 

from Eqs. (6) and (Al) for V{x) = V(y) with x G V, y = x G V. This invariant quantity Qp 
is the generalized version of Q in Eq. (2) for a general (smooth) similarity transform F, with 
alternative Qf defined in the same manner (replacing -0 by 0* in Eq. (A3)). 

In particular, we can express Eq. (A3) for basis functions <j>i(x), 02 (x) and 0f (y), 0f (y) of the 
solution spaces of Eqs. (6) and (Al), respectively, as 

1 


q F - = 

1 * ran 


2 i 


(p m (x)F\x)(p^(y) - 0' m (x)0;(i/) 


(AS) 
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in analogy with Eq. (8), with q^ ri (x, x) being spatially constant in T>. These generalized invariants 
now map the local basis (j>{x) in the original axis within a domain T> 3 x to the basis 4>f[x) hr 
the transformed axis within the image domain T> 3 x = y: 

4> f (x) = Q F <t>{x), (A6) 

where the mapping matrix Q F = Qf(%, x) has the same form as Q in Eq. (14) but with the q mn 
replaced by q ^ nM , as simply shown in the same manner as Eq. (13) from Eq. (12). 

In the present article we focus on the symmetry-induced mapping of a single basis between 
F-rnapped domains spanning the solution space of the Schrodinger equation in a fixed coordinate 
system. In other words, we demand that the basis functions 0f 2 in Eq. (A5) (or equivalently, 
the function ipp in Eq. (A3)) be solutions of the original equation (Eq. (6)), which is nontrivially 
the case only if Dp{y) = 0 f(u ) hi Eq. (Al). In other words, we here consider transformations 
which leave the Schrodinger equation invariant, so that the same basis <f> (or solution 0) at two 
points x,x can be used in Eq. (A5) (or in Eq. (A3)). Since this should hold for arbitrary locally 
F-symmetric potential (and thereby arbitrary 0 or 0), we demand that (cf. Eq. (A2)) 

F"(x) = 0, [F\x)} 1 2 3 4 5 6 7 8 9 = 1 => F(x) =ax + p, o = ±1, (A7) 

so that F is an isometry corresponding to the (local) inversion or translation transforms in 
Eq. (5). Further, local symmetry transforms of this type between finite domains T> and T> have 
a global limit with T> — M, which is not the case for an arbitrary (non-isometric) transformation 
F. The two particular cases of local inversion and translation symmetry (see Fig. 2 (b) and 
(c), respectively) entailed in Eq. (A7) then enable the recovery of stationary parity and Bloch 
eigenfunctions for the Hamiltonian in the global symmetry limit. 
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